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Abstract 

Starting from a model of an elastic medium, we derive equations of motion that are identical 
in form to Dirac's equation for a spin 1/2 particle with mass, coupled to electromagnetic and 
gravitational interactions. The mass and electromagnetic terms are not added by hand but emerge 
naturally from the formalism. A two dimensional version of this equation is derived by starting 
with a model in three dimensions and deriving equations for the dynamics of the lowest Fourier 
modes assuming one dimension to be periodic. Generalizations to higher dimensions are discussed. 
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I. INTRODUCTION 



Dirac's equation describes the behavior of particles with mass and spin and how they 
couple to the electromagnetic field. The usual form of Dirac's equation is 

(i^d^-m)m(x) = 

The electromagnetic field is introduced by the minimal coupling prescription^ <9 M — > D^, 
with 

where is the electromagnetic vector potential. Dirac's equation can be further coupled 
to gravity (at the classical level) using the prescription 3 

<9 M -> <9 M - r M 

and the equation then takes the for m 2 ' 3 ' 4 ' 5 ' 6 

- lF » - 4*1*0*) - ™*(*) = o (i) 

where is known as the spin connection, is the electromagnetic vector potential, and m 
is the mass. The gravitational coupling enters through the modified dirac matrices 7^ which 
satisfy the anticommutation relation 

{r,r}=i!T- 

and the operator 7 A '[c? M — T^} is (in the absence of electromagnetic interactions) the covariant 
derivative for spinor fields in a curved spaced. 

The above form of Dirac's equation describes the dynamics of the spinor field ^ when 
coupled to the scalar fields A^ and gravity. There are two additional equations which describe 
the dynamics of A^ and g^ u , these are the Einstein field equations 

Rfiv — -R = (2) 

and Maxwell's equations 

V^F"" = 4^*7"* (3) 
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The Equations (0) and (J3J) are collectively known as the Einstein-Dirac-Maxwell 
equations£&£. The subject of this paper is Equation (0). We will show that the equa- 
tions of motion of an elastic solid have the same form as Equation ((TJ with the mass end 
electromagnetic term emerging naturally from the formalism. 

II. ELASTICITY THEORY 

The theory of elasticity is usually concerned with the infinitesimal deformations of an 
elastic bod y 8 ^ 10 ' 11 ' 12 . We assume that the material points of a body are continuous and can 
be assigned a unique label a. For definiteness the elastic body can be taken to be a three 
dimensional object so each point of the body may be labeled with three coordinate numbers 
dj with i — 1, 2, 3. 

If this three dimensional elastic body is placed in a large ambient three dimensional space 
then the material coordinates ai can be described by their positions in the 3-D fixed space 
coordinates Xi with % = 1,2,3. In this description the material points di{x\, 22, ^3) are 
functions of x. A deformation of the elastic body results in infinitesimal displacements of 
these material points. If before deformation a material point a® is located at fixed space 
coordinates Xi,X2,£3 then after deformation it will be located at some other coordinate 
x[, x' 2 , x' 3 . The deformation of the medium is characterized at each point by the displacement 
vector 

/ 

which measures the displacement of each point in the body after deformation. 

It is the aim of this paper to take this model of an elastic medium and derive from it 
equations of motion that have the same form as Dirac's equation. 

We first consider the effect of a deformation on the measurement of distance. After our 
elastic body is deformed, the distances between its points changes as measured with the fixed 
space coordinates. If two points which are very close together are separated by a radius vector 
dxi before deformation, these same two points are separated by a vector dx\ = d.Xi + dui. 
The square distance between the points before deformation is then ds 2 = dx\ + dx\ + dx\. 
Since these coincide with the material points in the undeformed state, this can be written 
ds 2 = da\ + da 2 , + da 2 y The squared distance after deformation can be written^ ds 12 = 
dxf + dx' 2 + dx' 2 = J2i dxf = J2i(ddi + dui) 2 . The differential element dui can be written as 
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dui = J2i ada k , which gives for the distance between the points 



ds 
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da { + V ^dai 
i oai 



I 3u ■ du * du ■ du ■ 

E ddidcii + V —^-daida k + V —^dajdai + E E ^- L ^ 



ee(^+( 



E + 2( 4) ^Ma* 



da k dcii 



+ E 



dui dui 
daj dai. 



dakdai 



i 



where e' ik is 



, 1 / dui duk s-^ 
hk = 2 [da~k + ~da~ + ^ 



dui dui 
daj dai 



(4) 



; ^"^l ^^k j 

The quantity e' ik is known as the strain tensor. It is fundamental in the theory of elas- 
ticity. In most treatments of elasticity it is assumed that the displacements Uj as well as 
their derivatives are infinitesimal so the last term in Equation (jlj) is dropped. This is an 
approximation that we will not make in this derivation. 
The quantity 

dui dui 



9ik 



dui du k 

Oi,k + ~ h 



da k 

kk + 2 4 



da,; 



E 



da, da k 



(5) 



is the metric for our system and determines the distance between any two points. 

That this metric is simply the result of a coordinate transformation from the flat space 
metric can be seen by writing the metric in the form^ 



where 
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and J is the Jacobian of the transformation. Later in section fVTI we will show that the metric 
for the Fourier modes of our system is not a simple coordinate transformation. 
The inverse matrix (g lk ) = {gik)~ l is given by (g lh ) = (J _1 )(J _1 ) T where 



/ dai da\ dai \ 



J- 





dx' 2 


dx' 3 


da 2 


da 2 


da 2 


dx[ 


dx' 2 


dx' 3 


da 3 


da 3 


da 3 



\ dx'i d x 2 d x 3 / 



This yields for the inverse metric 



Jk 



Oik — 



dm 



dx k 



du k dui dui 
dxi , dxi dxk 



(6) 



(7) 



where is defined by 



1 / dui_ duk _ \- 
" lk ~ 2\dx k + d Xl ^ 



dui dui 
dxi dxi 



We see that the metric components involves derivatives of the displacement vector with 
respect to the internal coordinates and the inverse metric involves derivatives with respect 
to the fixed space coordinates. 



III. EQUATIONS OF MOTION 

In the following we will use the notation 

_ du^ 
U " v ~ dx v 

and therefore the inverse strain tensor is 

We will use the lagrangian method to derive the equations of motion for our system. Our 
model consists of an elastic solid embedded in a 3 dimensional euclidean space. 

In the following we work in the fixed space coordinates and take the strain energy as the 
lagrangian density of our system. This approach leads to the usual equations of equilibrium 
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in elasticity theoryS*^. The strain energy is quadratic in the strain tensor & v and can be 
written as 



iivap 



The quantities C^ ap are known as the elastic stiffness constants of the material. For an 
isotropic space most of the coefficients are zero and in 3 dimensions, the lagrangian density 
reduces to 



fc 12 ^ fc 13 ^ fc 33 



L = (A + 2/i) [e?! + e\ 2 + e§g] + 2A [e xl e^ + ene 33 + e 2 2e 33 ] + 4/x 

where A and /x are known as Lame constants^. 
The usual Lagrange equations, 

y. d ( dL \ 8L _ 
„ dx v \du pu J du p 

apply with each component of the displacement vector treated as an independent field vari- 
able. Since our Lagrangian contains no terms in the field u p , Lagrange's equations reduce 
to 

V — (-^-\ 

V dx v \du pv ) 

The quantity 

T J_fdL\ 



; , dxy \dtipv J 

is a vector and as such can always be written as the sum of the gradient of a scalar and the 
curl of a vector or 

V = V0 + V x A 
From this decomposition we can immediately conclude, 

V 2 = (9) 

We see therefore that the scalar quantity 4> in the medium obeys Laplace's equation. 



A. Physical Interpretation of <p 

To understand the physical origin of (f) we derive its form in the usual infinitesimal theory 
of elasticity. The advantage of the infinitesimal theory is that an explicit form of the vector 
V may be obtained. In the infinitesimal theory of elasticity the strain components Up U are 
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assumed to be small quantities and therefore the quadratic terms in the strain tensor are 
dropped and the strain tensor reduces to^ 



From this we see that the scalar in the infinitesimal theory is the divergence of the strain 
field a — V • u. It is an invariant with respect to change of coordinates and in general varies 
from point to point in the medium. This exercise exhibits the physical origin of which to 
lowest order in the strain components is the divergence of the strain field. 

In this work however will not make the infinitesimal approximation and we will work with 
the scalar and not a. In most of what follows, the exact form of is not important. It is 
only important that such a quantity exists and obeys Laplace's equation. 

B. Internal Coordinates 

The central results of this work will be given in sections IVl and fVT] where we will take one 
of our internal coordinates to be periodic and we will Fourier transform all quantities in that 
coordinate. We therefore need to translate the equations of motion V 2 = from the fixed 
space coordinates to the internal coordinates. For clarity, in the remainder of this text we 
change notation slightly and write the internal coordinates not as dj but as x\ and the fixed 
space coordinates will be unprimed and denoted 2j. Now using Ui = x\ — Xi we can write 




Using the above Lagrangian we obtain the explicit form 



where a = uu + u 2 2 + ^33 = V • u. 

Finally taking the divergence of V yields 



VV = 



d_ dx^j)_ 
dxi ~* dxi dx'j 




(10) 
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Equation (fTUj) relates derivatives in the fixed space coordinates X{ to derivatives in the 
material coordinates x\. As mentioned earlier, in the standard treatment of elastic solids 
the displacements itj as well as their derivatives are assumed to be infinitesimal and so the 
second term in Equation (jl(Jj) is dropped and there is no distinction made between the Xi 
and the x\ coordinates. In this paper we will keep the nonlinear terms in Equation (jl(Jj) when 
changing coordinates. Hence we will make a distinction between the two sets of coordinates 
and this will be pivotal in the derivations to follow. 

We will now demonstrate that Laplace's equation Q implies Dirac's equation. 



IV. CARTAN'S SPINORS 

The concept of Spinors was introduced by Eli Cartan in 1913^. In Cartan's original 
formulation spinors were motivated by studying isotropic vectors which are vectors of zero 
length. In three dimensions the equation of an isotropic vector is 

xl + x 2 2 + xl = (11) 

for complex quantities Xj. A closed form solution to this equation is realized as 

*i = £o-& %2 + and x 3 = -2^i ( 12 ) 

where the two quantities £j are 



. Xi — IX 2 . j. \—X\—lXi 

^° =± i—r~ and ei=± v — — ■ 

That the two component object £ = (£o, £1) is a spinor— can be seen by considering a rotation 
on the quantities v\ = x\— 1x2, and v 2 = —x\ — 1x2- If V{ is rotated by an angle a, 

Vi — ► Vi exp(za) 

then the spinor component £0 is rotated by a/2. It is clear that the spinor is not periodic 
in 2-7T but in 47r. A quantity of this type is a spinor and any equation of the form (jlljl has a 
spinor solution. 

Laplace's equation 

( d 2 d 2 d 2 \ 1 n 



dx\ dx\ dxi 
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* 9 ) 


d 




dx 2 J 
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dx 2 . 


l = - 

' dzi 



can be viewed as an isotropic vector in the following way. The components of the vector 
are the partial derivative operators d/dxi acting on the quantity <p. As long as the partial 
derivatives are restricted to acting on the scalar field it has a spinor solution given by 

= 2 ( ~ i~xZT ) = ~ ( 13 ) 

and 

^ = -\\ir + 'ir\ =ir ( 14 ) 

where 

z = X\ + ix2 and z\ = —x\ + ix 2 
and the "hat" notation indicates that the quantities £ are operators. The equations 

^° ~ dz 

and 

*_± 

are equations of fractional derivatives of order 1/2 denoted £o — D 1 ^ 2 and £i = D)J 2 ■ Frac- 
tional derivatives have the property thal^ 

£1/2^1/2 = d_ 

z z dz 

and solutions for these fractional derivatives can be written^ 

DU=^§-J\z-ty^{t)dt (15) 

The exact form for these fractional derivatives however, is not important here. The important 
thing to note is that a solution to Laplace's equation can be written in terms of spinors which 
are fractional derivatives. 

If we assume that the fractional derivatives |o and |i commute then we also have 

(loll) 2 = lololili 
d d 



dz dz\ 



1 ( d d \ ( d d 



4 \dxz dx 2 J \dxz dx 2 
4 y dx\ dx\ J 
IdxJ 
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Using this result combined with Equations ()13|) and (JTlJl we may write for the components 
of our vector 

-2£o£i (16) 



dxi 



dx 2 
and 



" >«„ 2 +l?) (17) 



" s a. 



dx 3 

This result gives the explicit solution of our vector quantities in terms of the spinor 
quantities 



A. Matrix Form 



It can be readily verified that our spinors satisfy the following equations 

> = 



and in matrix form 
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is equal to the dot product of the vector = d/dx^ with the Pauli spin matrices 



dx 3 dx 2 
d d 



(19) 
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are the Pauli matrices. 
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So Equation (fHJj) can be written 



E drft = 0. (20) 

where we have used the notation £ = £0. This equation has the form of Dirac's equation in 
3 dimensions. It describes a spin 1/2 particle of zero mass that is free of interactions. 

B. Relation to the Dirac Decomposition 

The fact that Laplace' equation and Dirac's equation are related is not new. However 
the decomposition used here is not the same as that used by Dirac. In the usual method, 
starting with the dirac equation (i'y^d^)^ = and operating with — ry M <9 M yields Laplace's 
equation for each component of the spinor field. In other words this method results in not 
one Laplace's equation but several (one for each component of the spinor). Conversely if one 
starts with Laplace's equation and tries to recover Dirac's equation one must start with 2 
independent scalars ( 4 in the usual 4 dimensional case) in order to derive the two component 
spinor equation (|2T)|) . 

What has been demonstrated in the preceding sections is that starting with only one 
scalar quantity satisfying Laplace's equation Dirac's equation for a two component spinor 
may be derived. Furthermore any medium (such as an elastic solid) that has a single scalar 
that satisfies Laplace's equation must have a spinor that satisfies Dirac's equation and such 
a derivation necessitates the use of fractional derivatives. 

The form of Equation (j2*U|) is relevant for a massless, non-interacting spin 1/2 particle. 
We will now demonstrate that if one of our internal coordinates is taken to be periodic a 
mass term as well as gravitational and electromagnetic interaction terms appear in Dirac's 
equation. 

V. TRANSFORMATION TO INTERNAL COORDINATES 

In section IVTl we will take the x' 3 coordinate to be periodic and we will derive equations 
for the Fourier components of our fields. Since the elastic solid is assumed to be periodic 
in the internal coordinates we need to translate our equations of motion from fixed space 
coordinates to internal coordinates. Using Equation (fTUj) we can rewrite Equation (j2"Uj) . as 
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1-1 

or 



(21) 



E = 



where d'^ = d/dx'^ and 7'^ is given by 



7** = 7 M + E (22) 
The anticommutator of these matrices is 

= {7', 7"} + EM7".7"} + E»«« (7°. 7"} + E"^"^{7",7' 9 } 

/3 a a/3 

a 

EE <T 

This shows that the gamma matrices have the form of the usual dirac matrices in a 
curved spaced To further develop the form of Equation (}2*T]) we have to transform the spinor 
properties of £. As currently written £ is a spinor with respect to the Xj coordinates not the 
x\ coordinates. To transform its spinor properties we use a similarity transformation and 
write £ = where S is a similarity transformation that takes our spinor in x M to a spinor 
in x' . We will not attempt to give an explicit form for S. We simply assume (similar to 
reference ^]) that this transformation can be effected by a real similarity transformation. 

We then have 

& l £ = Wz+s& ll i 

Equation (j2~Tj) then becomes 

Using (d'^S'^S = — S'~ 1 (9| t 5'). This can finally be written as 

W- r /Jf=° ( 23 ) 
12 



where T M = (d'^S^S and 7 M = S~ lr y'S. Equation has the form of the Einstein-Dirac 
equation in 3 dimensions for a free particle of zero mass. The quantity d 1 —Y^ is the covariant 
derivative for an object with spin in a curved spaced. In order to make this identification, 
the field r M must satisfy the additional equation^ 

^7 + - Y u r + = (24) 

where Tg u is the usual Christoffel symbol. We will now show that this equation does hold 
for this form of 17 



A. Spin Connection 

To show that Equation (J24j) holds we consider the equation <9„7 = where the vector 7 is 

3 

7 = E 7% 

and ejj. is a unit vector in the direction. Since 7 is a vector, then the quantity d u j = is 
a tensor equation. Therefore, in the primed coordinate system we can immediately write 

E + 7^) e; = 

/x=l 

where 7'^ = 7 M + J2 a §^"7° * s the expression of 7^ in the primed coordinate system. Using 
7 //x = S^S -1 , we have 

E (Wrs- 1 ) + (s^- 1 )^) e ; = 

or 

E (W^rs- 1 + sid'^ns- 1 + sr&s- 1 ) + (s^s- 1 )^) e ; = 0. 

Multiplying by S* -1 on the left and S on the right yields 

E (s-\d> u s)r + (KY) + w&s-^s + *fT%) e ; = 

Finally, using T M = (9 A1 5' _1 )5' and again noting that d^S^S = —S^d^S we have, 



We have just demonstrated that in the internal coordinates, the equations of motion of 
an elastic medium have the same form as the free-field Einstein Dirac equation for a massless 
particle in three dimensions. 
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B. Physical Content 

Thus far all of the transormations that have been obtained are "trivial" in the sense 
that they only result due to changing coordinates from the unprimed coordinates x^ to the 
primed coordinates x' . Changes of coordinates of course do not result in any new physical 
content. In particular the metric derived in Equation (J3J) does not lead to a curved space. 
The Riemann curvature tensor calculated from Equation (jSJ) is identically zero. Likewise the 
spin connection T M is due solely to a gauge transformation £ — > S^' and as such contains no 
physical content since it can be removed by transforming £' — ► S^ 1 ^. 

What we will demonstrate in the following sections is that for a system where one coor- 
dinate is periodic, the resulting 2 dimensional quantities are NOT trivial. In other words 
the metric that determines the dynamics of the Fourier components of £ does in fact lead 
to a curved space and the spin connection cannot be removed by a gauge transformation. 
Furthermore the introduction of the fourier components will generate extra terms in Equa- 
tion (|2"T|) that imply a series of equations relevant for particles with mass coupled to fields 
that can be associated with electromagnetism. We will show that in the low energy approxi- 
mation (ie a system in which only the lowest few modes are present) the equations of motion 
are identical in form to Equation (JIJ. 



VI. INTERACTING PARTICLES WITH MASS 

In this section we again consider a three dimensional elastic solid but we take the third 
internal dimension to be compact with the topology of a circle. All variables then become 
periodic functions of x' 3 and can be Fourier transformed. 

In preparation for Fourier Transforming we isolate the terms involving x' 3 and rewrite 
Equation (|2*T|) as, 




We first transform the partial derivatives of the u v in equation (|25j) to obtain 

[*■ k 

where u u/X) k is the k th Fourier mode of du u /dx^ and k = 2m /a with a the length of the 
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circle formed by the elastic solid in the x' 3 direction and i is an integer. Equation (|25|) now 
becomes, 



E 7 M f ^4,0 + E + "S^sl t 

fl=l \ v=l J 

+ 7 3 ^4,0 + E U v^kK + U33,kd^j £ 

Next we transform the spinor (noting that it is periodic in Ana), 
with g = 27cj/a and j an integer. This yields, 

Jx' 3 (k+q/2) 



ji T i 

? e 2 3 



EE e ' 



E 7 M fojw) + E «^<% + %/2w) e 9 /2 

+ 7 3 i(q/2)5 k>0 + E + i(q/2)u 33)k £ q/2 



u=l 



0. 



This equation is independently true for each distinct value of k + q/2 = m/2 or 2k + q = m 
with k, q, m an integer. Writing q = m — 2k yields finally, 



E 

k 



2_, 7 »«4,0 + 2^ <V,fc<t + * o U3 »> k S(m-2fc)/2 



l/=l 



+ 7 h 



3 / . (m — 2fc) 



4,o + E u »3,kd'„ + 



.(m - 2fc) 



U>33,k £(m-2k)/2 



(26) 







z i/=i z 

This is an infinite series of equations describing the dynamics of the fields £ m . This set of 
equations describes the dynamics of our elastic solid and contains the same information as 
Laplace's equation. 

So far no approximations have been made. In the next section we will demonstrate that 
if only the lowest modes are present, this reduces to an equation that is identical in form to 
Equation (JTJ). 



A. Spectrum of Lowest modes 

We now consider a theory in which only the lowest few modes in Equation (j2B| are present. 
We therefore keep only the m = 0, ±1/2 modes we obtain the following 3 equations, 



2/2 \ 2 

E 7 M k + E Uo + E ^3,o^ 7 3 £o = 

|U=1 V v=\ ) v=\ 



(27) 
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E 7 M I + E u ^,od' u + im 1/2 u 3fl)0 f 1/2 + l 3 im 1/2 {l + 1*33,0)6/2 

ji=l V v=l / 

2 2 

+ 7 3 E ^3,o%/2 + 7 3 E ^3,1^-1/2 = (2* 



u=l v=l 



E 7 M ( 9J, + E "^.o^, + wn_i /2 U3 M ,o £-1/2 + 7 3 ^-i/2(l + 1*33,0)^-1/2 

2 2 

+ 7 3 E ^3,0^-1/2 + 7 3 E ^3,-1^6/2 = (29) 

v=l v=l 

where = 2iri / a denotes the Fourier mode with i a half integer. These equations describe 
the dynamics of three fields £ an d the coupled fields £i/ 2 and £-1/2- The first equation 
(m = mode) describes the dynamics of a massless, free particle. We will not attempt to 
identify this mode with any physical particle but we simply note that in this approximation 
this equation is completely uncoupled from the m = ±1/2 modes and therefore its dynamics 
are independent and have no affect on these other modes. 

We now examine the equations describing £]y 2 and £-1/2 ■ These two equations can be 
combined by noting that for real fields, = u*_ h . The m = ±1/2 modes can now be 
combined into the single equation 

E 7 M + E Uvho&v + ^1/21*3^,0) * + 7 3 WH/2(1 + m 3 3,o)^ 

(1=1 \ v=l ) 

2 2 

±7 3 E ^3,0^* ±7 3 E ^3,1^** = 0, (30) 

v=l v=l 

where \l/ = £1/2 + £li/ 2 - To put this equation into a more recognizable form we multiply 
Equation (J3U)) by 7 3 from the left and define 

7' M = 7V±E7 3 7V/3,o- ( 31 ) 

(3=1 

These matrices with fi — 1,2 and v — 1, 2 satisfy the anticommutation relations 

2 

{V, l' 1 *} = <W + ( u ^,o + u ^,o) + E u ^fi u uPfi- 

13=1 

If we insist that our new matrices satisfy {7^,7"} = g^ v then we are led to define 

2 

9^ = V + (U m0 + U Vfi)0 ) + E u ^fi u u(3,0- (32) 

(3=1 
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This is the metric for our two dimensional subspace and it does not have the form of a simple 

coordinate transformation on a flat space metric like that of section [H] 

Equation (|50|) can now be rewritten as 

2 2 

E 7 /M + im 1/2 u 3flfi ) * - im 1/2 u 3fl) o 7 3 7 / V/3,o 1 I f + imi/ 2 (l + «33,o)£i/2 

M=l 0=1 

2 2 

+ E + E ^3,1^** = 0. (33) 

l/=l l/=l 

As we did for Equation (J25)l . in going from the a; to the x' coordinates, we assume that 

the spinor properties of £1/2 may be transformed using a real similarity transformation and 

writing £1/2 = S'£i/2- Transforming \I/ in this way and multiplying on the left by S^ 1 gives 

us the following form for the m = ±1/2 modes, 

2 2 

£ S-VS + zm 1/2M3 M,o + S-^S) * - m 1/2 u 3ft o E S'" V^o* (34) 

A*=l /3=1 

2 2 

+ zm 1/2 (l + 1/33,0)^ + E + S 1 - 1 ^)* + £ + S-^S)** 

We now examine each quantity in Equation (|34p. As before, we identify 7 M = S ^y'^S 
with the transformed gamma matrix and T M = (d' /1 S~ 1 )S with the spin connection. We also 
would like to identify the quantity = 777,11*3^0 in the first term with the electromagnetic 
potential and the third term in Equation (|3*4*j) as a mass term with m = 7774(1 + 1*33,0) which 
implies that the field 1*33,0 provides a mass for our \1/ particle. Let us further assume that the 
quantities are small compared to unity so that the second term may be neglected as being 
of order (ie we are now assuming that our medium undergoes only small deformations). 
If these identifications are made we can write Equation (|3*4"Jl in the final form 



]T r K + iT, -A,)*- (35) 



2 2 

+ 1 E - r„)* + * E - r„)**. 

f=l v=l 

Notice the formal similarity of this equation to Equation {1}. The first two terms have exactly 
the form of Dirac's equation for a spin 1/2 particle of mass m in curved space interacting with 
the electromagnetic vector potential A^. Note that the mass term and the electromagnetic 
potentials were not added by hand but emerged naturally from the formalism. The nature of 
the last two terms in Equation ()35|) are unknown. They don't appear in the usual statement 
of Dirac's equation and their implications are unknown. 
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Equation (|35p is the central result of this work. We have as yet not shown that the 
dynamics of the fields are consistent with their identification as the electromagnetic vector 
potential. To truly claim that the quantity is the electromagnetic potential it must be 
shown to satisfy Maxwell's equations. We believe however that the formal correspondence 
between Equation (J35j) and Dirac's equation is significant in its own right and will not, in this 
paper, pursue the question of whether Maxwell's equation or the Einstein Field equations 
are satisfied. Before concluding we note that although our derivation assumed that we were 
working in three dimensional space, the formalism extends to any number of dimension o 14 ! 17 . 
The major difference is that in the three dimensional case, we were able to find an explicit 
solution for the components of a spinor in terms the components of the vector <9 M . An explicit 
solution might not exist in general. Nevertheless it can be shown 1 ^ that the quadratic form in 
Laplace's equation implies the existence of a multicomponent spinor £ satisfying a dirac-like 
equation in any dimension. 



VII. CONCLUSIONS 



We have taken a model of an elastic medium and derived an equation of motion that 
has the same form as Dirac's equation in the presence of electromagnetism and gravity. 
We derived our equation by using the formalism of Cartan to reduce the quadratic form 
of Laplace's equation to the linear form of Dirac's equation. We further assumed that one 
coordinate was compact and upon Fourier transforming this coordinate we obtained, in a 
natural way, a mass term and an electromagnetic interaction term in the equations of motion. 
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